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ABSTRACT 

An  EPR  study  of  copper  halide  tetrazole  complexes  and  x-irradiated 
strontium  acetate  hemihydrate  has  been  made.   The  powder  spectra  of  some 
of  the  tetrazole  complexes  exhibited  temperature  dependent  narrowing  and 
broadening.   The  one-electron  population  of  the  copper  d  orbital  involved 
in  bonding  to  2-methyl-5-amino  tetrazole  has  been  found  to  be  0.56  and 
0.47,  respectively  for  the  chloride  and  bromide  complexes.   The  g  value 
for  these  two  complexes  are  consistent  with  the  proposed  planar  poly- 
meric structure  for  2-substituted  tetrazole-copper  complexes.   The 
species  present  in  single  crystals  of  strontium  acetate  hemihydrate  ir- 
radiated and  observed  at  liquid  nitrogen  temperature  to  -100  C  has  a 

CH 
spectrum  of  four  lines  of  1:3:3:1  intensity  ratio  and  a   3  =  11.0  gauss. 

Temperature  dependent  exchange  behavior  of  the  spectrum  for  this  species 

is  observed  and  attributed  to  a  rotating  methyl  group  which  was  found  to 

have  AlT  =  3.0  +  0.3  -—■     and  AS^  =  6.0  +  0.3  cal/mole  deg  for  internal 
—     mole  — 

rotation.   The  species  present  at  room  temperature  has  an  eight  line 

CH 
spectrum  consisting  of  two  1:3:3:1  intensity  ratio  patterns  with  a   3  = 

25.2  and  a   =  21.3  gauss.   Molecular  orbital  calculations  using  the  CNDO 

approximation  were  performed  and  possible  fragments  responsible  for 

these  spectra  are  suggested. 
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I.   INTRODUCTION 

The  work  reported  in  this  thesis  consists  of  two  parts.   The  first 
part,  Sections  II  through  VI  inclusive,  contains  the  results  of  research 
on  copper  halide- tetrazole  complexes.   The  second  part,  which  encompasses 
the  rest  of  the  thesis,  examines  the  problem  of  the  paramagnetic  frag- 
ments produced  in  strontium  acetate  hemihydrate  by  x-irradiation. 

The  tetrazole  work  was  suggested  by  Dr.  H.  B.  Jonassen  of  the 
Chemistry  Department  of  Tulane  University,  New  Orleans,  Louisiana,  to 
provide  information  to  supplement  that  already  gathered  by  other  means 
to  study  the  structure  of  the  copper  halide  tetrazole  complexes  and  the 
coordination  of  the  tetrazoles  to  the  copper  (II)  in  the  complex.   The 
radiation  damage  work  on  strontium  acetate  was  begun  at  the  conclusion 
of  the  copper  tetrazole  work  with  the  purpose  of  obtaining  information 
about  the  mechanism  of  radiation  damage  in  crystals.   This  work  has 
resulted  in  the  discovery  of  free  radical  fragments  not  previously  re- 
ported for  a  crystal  of  this  type  [25]. 

For  each  section,  the  thesis  is  written  so  that  the  reader  is  given 
the  theory  and  background  necessary  to  understand  and  appreciate  the 
reported  results  in  the  context  of  current  research  in  the  areas  of 
transition  metal-tetrazole  complexes  and  radiation  damage  studies  of 
organic  crystals. 


II.   ESR  OF  TRANSITION  METALS 

Electron  spin  resonance  (ESR)  measures  the  absorption  spectra  which 
arise  from  transitions  between  energy  states.   The  energy  states  are  the 
result  of  an  interaction  of  the  unperturbed  energy  states  of  the  para- 
magnetic species  with  a  magnetic  field.   In  this  section  the  theory  of 
these  states  and  their  description  by  a  spin  Hamiltonian  will  be  dis- 
cussed.  The  theory  of  the  g  tensor  for  transition  metals  as  developed 
by  Carrington  and  McLachlan  [l]  will  also  be  considered. 

Elementary  quantum  mechanics  provides  the  following  relationship 
for  calculating  the  average  value  of  a  property  of  a  quantum  mechanical 
system: 

<q>  =  ^  f'e^dT 

where  9  is  an  operator  for  the  property  Q  whose  average  value  is  to  be 
calculated  and  f  is  the  spatial  wave  function  for  the  system.   ESR  is 
used  to  determine  the  energy  levels  of  a  system.   An  extension  of  the 
above  relationship  to  a  system  with  several  energy  levels  can  be  made 
without  too  much  difficulty.   The  operator  for  the  energy  is  called  the 

spin  Hamiltonian  which  is,  for  one  electron: 

— ♦  -♦ 

K     =PH  •  g  •"s-Sv.^H  •  1.  +E"S  -A.  *  I. 

S  I     '  1         1  11 

In  essence,  the  terms  in  the  Hamiltonian  describe  the  contribution 
to  the  total  energy  which  are  made  by  interactions  of  the  electrons  and 
nuclei  with  the  external  magnetic  field  and  with  one  another.   The  first 
term  is  the  electron  Zeeman  term.   It  represents  the  Zeeman  energy,  the 
energy  of  interaction  of  the  electron  spin  moment  with  the  external 
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magnetic  field.   |3  is  the  Bohr  magneton,  S  is  the  electron  spin  and  g 
is  the  spectroscopic  splitting  tensor  which  may  well  be  anisotropic. 
The  second  term  is  a  summation  of  the  Zeeman  energies  of  the  nuclei. 
I.  is  the  spin  of  the  ith  nucleus  and  y.    is  the  gyromagnetic  ratio  of 
the  ith  nucleus.   The  third  term  represents  the  summation  of  the  energy 

of  the  hyperfine  interaction  between  the  electron  spin  and  the  spin  of 

— * 
the  ith  nucleus.   A.  is  the  hyperfine  interaction  tensor.   In  this 

section  only  the  first  term  in  the  Hamiltonian  will  be  considered: 

K     =  6H  •  g  •  "S  =  (H  H  H  )   /g    g   g 
s  x  y  z    /  xx  xy  xz 

egg 

yx  yy  yz 
8zx  8zy  Szzj 

In  the  absence  of  spin-orbit  coupling  the  g  tensor  is  isotropic  and 
S  represents  the  true  spin  of  a  free  electron.   However,  consider  the 
case  where  g  is  anisotropic,  since  this  is  the  case  which  is  almost  al- 
ways encountered  with  transition  metal  complexes.   It  naturally  follows 
from  above  that  an  anisotropic  g  tensor  implies  that  S  cannot  possibly 

represent  the  true  spin.   Instead,  for  anisotropic  g, 

-♦   -» 

K     =   0  H  •  g  *  S 
s 

A 

where  S  is  called  the  fictitious  spin. 

The  reader  might  recall  the  fact  that  in  the  presence  of  a  magnetic 
field,  there  is  an  interaction  between  H  and  L,  the  orbital  angular 
momentum  of  the  electron.   In  this  case,  the  true  form  of  the  Hamiltonian 
for  the  interaction  of  an  unpaired  electron  with  a  magnetic  field  is 


K      -  0H  •  L  +  g  3H  •  S 
s  e 

where  e   is  the  free  electron  g  =  2.0023.   The  two  forms  of  H  are 
°e  s 
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equivalent.  This  fact  is  used  to  develop  some  mathematical  relation- 
ships which  will  help  in  the  understanding  of  the  anisotropy  of  the  g 
tensor. 

It  can  be  shown  that  the  only  way  an  odd  electron  can  acquire  or- 
bital angular  momentum  is  through  the  effect  of  spin-orbit  coupling, 
which  is  represented  by 

-♦   -♦ 
C  L  •  S 

where  C,     is  called  the  spin-orbit  coupling  constant,  the  value  of  which 

may  be  extracted  from  atomic  spectra. 

First  order  perturbation  theory  allows  the  calculation  of  modified 

— ♦   — » 
ground  state  wave  functions  <  +  1  and  <  -  1.    £  L  •  S  is  the  pertur- 
bation operator  which  mixes  the  ground  state  wave  functions  (\|j  fa\    and 
(\J;  ,3 1  with  excited  states.  \<y)   and  |(3)  are  the  spin  functions  (for 
notation,  and  details  of  development,  see  [l]  ). 

The  fictitious  spin  can  now  be  defined  as  follows: 

*  i  +">  =  %i  +">       s  i  +>  =  %i  -  > 

Z  A. 

s  I  -  )  =  -%|  -  )      s J  +  )  =  %i|  -  > 

&  y 

In  a  magnetic  field  along  the  i  axis  (i  =  x,y,z) 

K     -  aH.  -(g.  S  +  g.  S  +  g.  S  )  . 
s      l    ix  x    ly  y    iz  z 

Setting  up  the  matrix  representation  of  this  form  of  3C  in  the  basis 

+  ),    -)  and  comparing  it  with  the  matrix  representation  of 

K     =   (3H  *  L  +  g  f3H  •  S  in  the  same  basis,  it  is  found  that: 
s  °e 
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e       -  g     -   2       E 


<0|L    |n><n|L    |  0> 


x 


'xx   °e       n       AE 

n 


(OIL  |n><n  L  I  0> 

g     -g   -  2    s       r„      r 


yy   e     n     ae 


n 


g.  -  g  -  2   S 


<0|L    |n>(n|L    I  0> 


3zz        °e  „  AE 

11  n 

where  AE   is  the  difference  in  energies  between  the  ground  state  elec- 
tronic level  and  the  nth  electronic  level.   The  value  of  AE  may  be 

n 

extracted  from  optical  spectra. 
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III.   EXPERIMENTAL  SPECTRA  AND  THEIR  ANALYSIS 

A.   SINGLE  CRYSTAL  METHOD 

Of  the  methods  available  for  studying  paramagnetic  species  in  solids, 
the  single  crystal  method  provides  the  best  results.   The  species  to  be 
studied  is  regularly  oriented  relative  to  the  crystal  axes,  and  this 
provides  the  opportunity  to  rotate  the  crystal  in  the  magnetic  field  and 
measure  the  anisotropic  interactions,  giving  valuable  information  about 
structure  of  the  paramagnetic  species. 

The  type  of  spectrum  obtained  is  dependent  upon  the  orientation  of 
the  external  magnetic  field  with  respect  to  the  paramagnetic  species, 
and  the  number  and  type  of  different  species  in  the  crystal.   In  addi- 
tion, there  may  be  more  than  one  magnetically  distinct  site  for  a  parti- 
cular species  in  the  lattice.   Fortunately,  there  is  usually  a  parti- 
cular orientation  for  the  magnetic  field  which  results  in  a  spectrum 
which  is  relatively  easy  to  interpret  and  allows  the  experimentalist 
to  extract  the  desired  information.   This  information  enables  the  experi- 
mentalist to  find  the  orientation  of  the  principal  axis  system  of  the 
species  with  respect  to  the  crystal  axis  system.   It  is  in  the  principal 
axis  system  that  the  A  and  g  tensors  are  diagonal.   Since  each  tensor 
has  three  diagonal  components,  there  are  three  orthogonal  directions  for 

the  magnetic  field  along  which  the  values  of  g   ,  g   ,  and  g   may  be 
°  °  °xx   yy       zz 

extracted.   A   ,  A   ,  and  A   may  sometimes  also  be  extracted, 
xx   yy       zz 


14 


B.  THE  POWDER  AND  FROZEN  SOLUTION  METHODS 

When  single  crystals  of  the  material  to  be  studied  cannot  be  grown, 
or  when  the  information  cannot  be  extracted  from  single  crystal  studies, 
then  the  powder  or  the  frozen  solution  "glass"  method  is  utilized.   In 
both  the  powder  and  glass,  the  paramagnetic  species  are  randomly  orient- 
ed with  respect  to  the  magnetic  field  and  absorb  over  a  wide  range  of 
magnetic  field.   The  glass  provides  a  method  for  diluting  the  species 
in  order  to  reduce  the  line  broadening  effects  which  are  caused  by  the 
interaction  of  neighboring  complexes. 

C.  LINE  WIDTHS 

There  are  three  factors  which  McGarvey  [2~|  considers  that  can 
contribute  to  the  width  of  ESR  absorption  lines.   They  are: 

(a)  Anisotropy  of  g  factors  and  hyperfine  interactions  or  spin- 
interactions  . 

(b)  Interaction  with  magnetic  dipoles  of  neighboring  electronic 
and  nuclear  spins. 

(c)  Spin-lattice  relaxation  times. 

The  first  factor  is  important  for  powder  samples  in  which  the  small 
crystals  are  randomly  oriented.   In  this  case,  the  method  needed  to 
obtain  narrower  lines  is  to  grow  single  crystals.   If  this  cannot  be 
done,  then  it  is  sometimes  possible  to  extract  the  information  from  the 
shape  of  the  broad  powder  spectrum. 

The  second  factor  can  often  be  a  significant  one,  since  magnetic 
fields  produced  by  dipoles  over  distances  of  a  few  angstroms  are  very 
large.   These  large  internal  magnetic  fields  coupled  with  the  external 
magnetic  field  may  cause  the  local  field  at  the  site  of  the  paramagnetic 
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species  to  vary  greatly.   Since  the  field  of  a  magnetic  dipole  is  in- 
versely proportional  to  the  third  power  of  the  distance,  the  best  way  to 
reduce  this  broadening  effect  is  to  place  the  species  in  an  isomorphic 
diamagnetic  crystal  or  to  use  the  glass  method. 

The  length  of  time  associated  with  the  relaxation  of  an  excited  spin 
state  to  a  thermal  equilibrium  spin  state  given  by  Boltzmann  statistics 
is  called  the  spin-lattice  relaxation  time.   The  relaxation  is  accom- 
plished by  a  transfer  of  energy  between  the  magnetic  spins  and  the 
vibrational  energy  of  the  crystal  lattice  to  which  the  spins  are  coupled. 
For  systems  where  the  spins  are  strongly  coupled  to  the  vibrational 
modes,  the  lifetime  of  a  given  magnetic  state  is  short,  resulting  in 
an  uncertainty  in  the  energy  which  manifests  itself  as  a  broad  absorp- 
tion line.   This  broadening  mechanism  is  strongly  dependent  on  the 
temperature,  so  that  the  lines  broadened  in  this  manner  can  be  sharp- 
ened by  lowering  the  temperature  of  the  sample.   Extremely  short  relax- 
ation times  often  occur  when  the  paramagnetic  species  has  an  electronic 
excited  state  just  above  the  ground  state.   In  cases  of  this  sort,  it 
is  often  necessary  to  go  to  liquid  helium  temperature. 
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IV.   BACKGROUND  ON  TRANSITION  METAL- TETRAZOLE  COMPLEXES 

Tetrazoles  are  f ive-membered,  heterocyclic  ring  compounds  which 
contain  one  carbon  and  four  nitrogen  atoms  linked  by  three  single  and 
two  double  bonds.   The  parent  compound  may  exist  in  tautomeric  forms  I 
and  II.   It  was  reported  by  Lounsburg  [3]  that  977o  of  tetrazole  is  of 
the  form  I : 

(1)N C(5) 


N-      C 


(2)N_       N(4)  H—  ^ 

VN: 
(3) 

I  II 

The  tetrazole  ring  is  unusual  in  that  the  only  two  points  of  sub- 
stitution are  in  positions  1  or  2  and  5.   Chemical  investigations  have 
shown  that  the  tetrazoles  are  nucleophilic  reagents  and  that  the  nucleo- 
philic  character  varies  with  the  nature  and  position  of  the  groups  sub- 
stituted on  the  ring. 

Jarvis  [4]  investigated  the  crystal  structure  of  a  1:1  complex  of 
copper  (II)  chloride  and  1,2,4  triazole.   He  reported  the  unit  cell 
dimension  to  be 

a  =  6.81,    b  =  11.39,    c  =  7.13  A  ;    0  =  96°  58' 
The  space  group  is  I2/c  and  contains  four  units  of  C-N-H  *  CuCl  .   The 
copper  atom  has  a  distorted  octahedral  coordination  group  consisting  of 
two  N  at  1.98  A,  two  CI  at  2.34  A,  and  two  CI  at  2.77  A.   The  structural 
unit  is  an  infinite  chain  in  which  octahedral  groups,  joined  by  sharing 
edges,  are  also  linked  by  triazole  molecules. 
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"N- 


I  I 

Nv      N 

V 

H 
1,2,4  TRIAZOLE 

Jonassen  [5]  reported  the  preparation  of  a  number  of  copper  halide 
tetrazole  complexes  of  differing  compositions.   He  determined  that  the 
copper  halide  complexes  of  1-substituted  tetrazoles  have  the  composi- 
tion Cu  X_T  and  postulated  a  structure  analagous  to  the  triazole  struc- 
tural unit  for  them.   The  copper  halide  complexes  of  2-substituted 
tetrazoles  have  the  composition  CuX  T.   The  proposed  structures  of  the 
1  and  2-substituted  tetrazole  complexes  of  copper  (II)  halides,  respec- 
tively,  are  shown  in  the  top  and  bottom  of  Fig.  1.   Jonassen  also 
carried  out  a  reduction  reaction  on  several  of  the  copper  tetrazole 
complexes  with  acetone  and  discovered  that  the  product  of  the  reaction 
was  always  a  CuXT,  regardless  of  the  composition  of  the  starting  material 
He  reasoned  that  in  CuXT,  polymeric  CuCl  chains  are  present  with  bridg- 
ing tetrazole  groups  forming  four  coordinated  copper  (I)  species. 

Garber  and  Brubaker  [6],  reporting  the  results  of  their  work  on 
bis (1-methy 1-5-tetrazoy 1)  and  bis (1-cyclohexy 1-5- tetrazoy 1)  nickel  (II), 
stated  that  the  reflectance  spectra  of  these  complexes  indicated  that 
the  nickel  is  in  an  octahedral  environment.   They  also  concluded  that 
the  insolubility  of  these  complexes  in  all  common  solvents  suggest  a 
polymeric  structure. 
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Brubaker  [7]  carried  out  continuous  variation  experiments  on  5- 
aminotetrazolato-copper  (II)  and  showed  that  there  is  a  1:2  interaction 
of  copper  (II)  with  the  tetrazole  according  to  the  mechanism: 

HT  =  H  +  T"  K  =  6.8  x  10~7 

a 

Cu2+  +  2T~  =  CuT 

He  determined  the  formation  constant  of  the  association  reaction  to  be 

12 
K  =  1.3  x  10   .   Infrared  analysis  of  this  complex  indicated  that  co- 
ordination of  the  tetrazole  to  the  copper  (II)  does  not  involve  the 
amino  group.   Jonassen  [5],  reporting  on  this  same  complex,  stated  that 
the  tetrazole  nitrogens  form  dative  bonds  with  the  copper  (II)  above  and 
below  the  plane,  thus  forming  a  bridged  structure. 

Bowers  and  Popov  [8]  worked  on  pentamethy lene  tetrazole  complexes 
of  copper  halide  and  compared  the  results  with  work  carried  out  pre- 
viously on  Cu(PMT) ,(C10. )„.   They  found  that  PMT  is  not  a  strong  donor 

b     4  2 

ligand  and  that  at  most,  two  PMT  coordinate  with  Cu(II)  in  the  halide 
complex  as  compared  with  six  PMT  with  the  perchlorate  complex.   In  a 
few  cases,  they  discovered  solubilization  of  complexes  when  the  solvent 
had  strong  donor  properties.   They  found  an  abnormally  low  magnetic 
moment  for  the  complex  which  is  indicative  of  copper-copper  bonding, 
but  they  concluded  that  their  data  did  not  allow  an  unambiguous  identi- 
fication of  such  bonds.   They  also  concluded,  from  magnetic  and  spectral 
evidence,  that  the  copper  (II)  in  Cu(PMT)X   is  in  an  octahedral  environ- 
ment . 

Kuska,  D'ltri  and  Popov  [9]  investigated  the  powder  spectrum  of 

Cu(PMT) , (CIO. )0  and  from  the  observable  nitrogen  superhyperf ine  splitting 
b    4  z 

and  comparison  with  optical  and  ESR  data  for  other  copper  complexes, 
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Y-N-C 


FIGURE  1.   Proposed  structures  of  1-substituted  tetrazoles  (top)  and 
2-substituted  tetrazoles  (bottom)  after  Jonassen  [5], 
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concluded  that  the  bonding  of  the  PMT  ligand  is  80%  ionic.   They  also 

concluded  that  the  PMT  nitrogen-copper  (II)  orbitals  have  a  greater 

2 

percentage  of  s  character  than  the  assigned  sp  hybridization. 
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V.   MATERIALS  AND  EQUIPMENT 

A.   OBTAINING  THE  DATA 

The  copper  tetrazole  complexes  studied  were  supplied  by  Dr.  R.  A. 
Henry  of  the  Chemistry  Division,  Naval  Weapons  Center,  China  Lake, 
California.   Their  purity  was  checked  by  Galbraith  Laboratories,  Inc., 
Knoxville,  Tennessee.   The  analytical  data  indicated  that  the  samples 
were  of  sufficiently  high  purity  for  this  study  [5]. 

The  ESR  spectra  were  obtained  with  a  Varian  V-4502-06  X-band  ESR 
spectrometer  utilizing  100  kc/sec  modulation  with  a  nine  inch  magnet 
and  field-dial  regulation.   The  cavity  was  a  conventional  reflection 
type  cavity,  equipped  with  a  Varian  variable  temperature  control  unit 
which  allowed  temperature  regulation  down  to  90  K.   For  liquid  helium 
work,  the  same  type  of  spectometer  modified  to  accomodate  a  liquid 
helium  dewar  and  pumping  system  was  used.   This  system  allowed  obser- 
vation of  the  sample  at  3  K. 

The  powder  spectra  were  obtained  by  placing  the  powder  in  4  mm 
O.D.,  3  mm  I.D.  quartz  tubing  sealed  at  one  end.   The  frozen  solution 
spectra  work  was  accomplished  in  a  similar  manner  by  filling  the  sealed  " 
quartz  tubing  with  a  solution  of  the  tetrazole  complex  and  then  freez- 
ing it  with  the  temperature  control  unit.   The  single  crystal  work  was 
accomplished  by  rotating  the  sample  on  the  gear  wheel  of  a  goniometer, 
an  instrument  which  gives  the  experimenter  the  freedom  to  orient  the 
crystal  about  two  axes,  the  axes  vertical  to  and  horizontal  with  the 
pole  faces  of  the  magnet. 


22 


The  data  was  taken  from  the  EPR  spectra  for  the  powder  and  frozen 
solution  method  by  using  the  baseline  of  the  spectrum  as  a  reference 
and  measuring  the  linewidths  and  coupling  constants  with  a  ruler  gradu- 
ated in  hundredths  of  an  inch.   The  data  for  the  single  crystal  work  was 
taken  from  the  oscilloscope  display  of  the  absorption  spectrum.   The 
orientation  of  the  single  crystal  was  changed  with  the  goniometer  while 
observing  the  scope  display,  and  data  taken  when  the  absorption  line 
reached  an  extreme  or  inflection  point  of  the  magnetic  field  at  which 
resonance  was  occurring.   The  resonance  fields  were  read  off  the  field 
dial  by  centering  the  absorption  line  on  the  crosshairs  of  the  oscillo- 
scope display. 

The  calibration  of  the  field  dial  with  the  linear  scale  of  the  re- 
corder was  accomplished  by  inserting  a  solid  mixture  of  0.1%  MnCl  in 
MgO  into  the  cavity  and  comparing  the  field  dial  readings  and  linear 
measurements  of  the  coupling  constant  from  the  recorder  with  the  coupl- 
ing constant  as  determined  with  an  NMR  probe. 

For  the  single  crystal  work,  each  piece  of  data  consisted  of  two 
angles  of  rotation,  a  resonance  field,  and  a  linewidth.   The  data  was 
cataloged  with  a  stereographic  net,  and  the  three  orientations  at  which 
data  was  taken  for  each  single  crystal  were  confirmed  to  be  orthogonal 
to  one  another  on  the  stereographic  net. 

B.   CATALOGING  THE  DATA 

Suppose  a  paramagnetic  single  crystal  is  to  be  studied.   This  crystal 
has  a  definite  form  which  allows  for  definition  of  a  three  dimensional 
orthogonal  crystal  axis  system  relative  to  the  crystal  faces.   The 
experimentalist  has  the  job  of  finding  the  orientation  of  the  principal 
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axis  system  with  respect  to  the  crystal  axis  system.   Consider  this  pro- 
cess from  a  different  viewpoint.   Reorienting  the  crystal  in  a  fixed 
magnetic  field  is  equivalent  to  reorienting  the  magnetic  field  with 
respect  to  a  fixed  crystal  axis  system. 

Consider  the  crystal  axis  system  X  Y  Z  in  Fig.  2.   It  is  easy  to  see 
that  all  points  in  space  may  be  convered   by  rotating  a  given  vector 
through  all  angles  of  6  and  (j>.      Suppose  the  experimenter  wants  to 
collect  data  from  the  paramagnetic  crystal  over  several  magnetic  field 
orientations.   The  stereographic  net  provides  a  convenient  method  of 
cataloging  the  information  obtained.   Each  value  of  9  and  <j>   corresponds 
to  one  point  on  the  stereographic  net.   Figure  2  gives  the  three  dimen- 
sional picture  of  a  principal  axis  system  X'  Y'  Z'  relative  to  the 
crystal  axis  system  X  Y  Z  .   The  representation  of  the  same  points  on 
the  stereographic  net  appears  directly  above  for  comparison.   The 
stereographic  net  covers  all  points  over  a  hemisphere.   Because  the 
ESR  spectrum  remains  unchanged  by  a  reversal  of  the  magnetic  field 
direction,  corresponding  points  on  the  two  hemispheres  are  equivalent. 

The  procedure  for  finding  X'  Y1  Z'  consists  of  taking  data  over 
several  orientations  of  the  magnetic  field  and  cataloging  the  data  with 
the  stereographic  net.   Once  this  is  done,  a  reasonable  guess  is  made 
at  the  three  orthogonal  points  on  the  net  which,  in  consideration  of 
the  data  distribution  on  the  net,  seem  most  likely  to  represent  X1,  Y1, 
and  Z'.   The  initial  guesses  are  then  fed  into  a  computer  program  for 
further  refinement. 

The  stereographic  net  thus  provides  a  convenient  method  for  catalog- 
ing data  and  providing  the  experimenter  with  an  overall  picture  of  the 
data  distribution  which  leads  to  an  initial  guess  at  the  orientation  of 
the  principal  axis  system  relative  to  the  crystal  axis  system. 

24 


Wolff   Net 


FIGURE     2 


z 


VI.   EXPERIMENTAL  WORK  AND  RESULTS 

A.   FROZEN  SOLUTION  SPECTRA 

Frozen  solution  work  was  undertaken  in  order  to  study  those  com- 
plexes which  did  not  grow  single  crystals  or  for  which  the  powder  spectra 
were  unsatisfactory.   The  following  organic  solvents  were  tried: 

Isopentane 

Glycerine 

Benzene 

Chloroform 

Diethyl  ether        and  mixtures  of  these 

Methyl  ethyl  ketone 

Acetone 

Ethanol 

Dimethylformamide 

A  small  amount  of  the  powder  was  placed  in  a  beaker  and  the  sol- 
vent added  at  room  temperature.   Upon  very  careful  attempts  to  dis- 
solve the  complexes  in  the  solvents,  it  was  discovered  that  in  almost 
all  attempts  either  the  complex  did  not  dissolve  at  all  or  what  at 
first  appearances  seemed  to  be  dissolution  of  the  complex  in  the  sol- 
vent was  in  actuality  only  a  fine  dispersion  of  the  powder  in  the 
solvent.   A  few  of  the  complexes  did  dissolve  to  a  very  slight  degree, 
however,  the  spectra  which  were  obtained  in  the  glass  at  113   K  were 
impossible  to  interpret  due  to  the  small  signal  to  noise  ratio  and 
other  complicating  effects. 

The  only  solvent  which  did  dissolve  all  of  the  copper  halide  tetra- 
zole  complexes  was  dimethylformamide.   The  copper  bromide  tetrazoles 
gave  quite  simple  spectra  in  DMF  whereas  the  spectra  of  the  copper 
chloride  tetrazoles  showed  several  overlapping  peaks  which  was  indicative 
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of  more  than  one  distinct  position  in  the  frozen  lattice  of  dimenthylfor' 
mamide. 

Three  typical  spectra  of  frozen  DMF  solutions  of  copper  bromide 
tetrazoles  are  shown  in  Fig.  3.   The  top  spectrum  is  that  of  copper  bro- 
mide dissolved  in  DMF.   It  was  used  as  a  reference  with  which  to  compare 
the  other  spectra.   The  other  two  spectra  exemplify,  respectively,  the 
phenomena  of  broadening  of  the  spectral  lines  and  shifting  of  the 
spectral  lines.   The  middle  spectrum  is  l-benzyl-5-amino  tetrazole  with 
copper  (II)  bromide  in  DMF,  and  the  bottom  spectrum  is  1,5  tetramethy- 
lene  tetrazole  with  copper  (II)  bromide  in  DMF. 

The  fact  that  the  complexes  did  not  dissolve  to  an  appreciable 
extent  in  all  but  DMF  substantiates  the  conclusion  of  previous  workers 
that  the  complexes  have  a  polymeric  structure  which  does  not  break  down 
easily.   The  dissolution  in  DMF  seems  to  imply  that  DMF  coordinates 
more  strongly  to  copper  than  do  the  tetrazoles.   The  data  in  Table  1 
leads  to  a  partial  confirmation  of  this  fact,  and  also  allows  for  some 
other  conclusions. 

Previous  workers  [9]  concluded  that  the  bonding  of  the  PMT  ligand 
to  copper  is  80%  ionic  and  that  PMT  is  a  poor  donor.   Data  in  Table  1 
indicates  that  the  PMT  ligand  has  been  replaced  by  DMF  and  that  the 
same  thing  has  happened  for  some  of  the  other  complexes.   This  is  not, 
however,  the  case  for  all  the  complexes  as  can  be  seen.   For  the  two 
complexes  which  exhibit  values  of  AH  and  A  different  from  those  listed 
for  the  copper  bromide  case,  it  must  be  concluded  that  the  tetrazole 
ligands  have  not  been  completely  replaced  by  the  DMF,  but  rather,  a 
type  of  coordination  still  exists  between  the  copper  and  the  tetrazoles. 
The  difference  in  AH  values  is  due  to  the  difference  in  the  ligand 
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FIGURE     3 
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COMPLEX 


AH  (gauss )_ 


AVERAGE  A  (gauss)* 


CuBr2 

2  phenyl 

1  Me  thy 1-5 -ethyl 
1,5  pentamethylene 

2  methyl-5-amino 
1  benzy 1-5-amino 
1  methyl-5-amino 
1,5  tetramethylene 
1,5  trimethylene 


448  +6  119  +  6~ 

442  +6  119  +  6" 

448  +6  116  +  6' 

448  +8  119  +  8' 

lines  too  broad  to  be  observed 

lines  too  broad  to  be  observed 

lines  too  broad  to  be  observed 

416  +  8  139  +  8' 

412  +8  138  +  8' 


'^Superscript  indicates  the  number  of  A  values  which  were  averaged. 

In  cases  where  less  than  three  values  were  averaged,  the  remaining 
A  values  were  lost  to  line  broadening  effects. 


-*  A  *- 


TABLE  1.   CuBr  -Tetrazole  Complexes  in  Frozen  Solution  of  Dimethylfor- 
mamiae 
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splitting  of  the  electronic  levels  of  the  copper  ion.   Since  the  AH 
values  are  inversely  proportional  to  the  ligand  splitting,  AE  (see 
Introduction),  larger  AH  values  imply  a  smaller  splitting  by  the  ligand, 
The  A  values  are  a  measure  of  the  localization  of  the  electron  in  the 
copper  orbital.   The  larger  the  A  value,  the  more  the  electron  is 
localized. 

B.   POWDER  SPECTRA 

Due  to  the  lack  of  success  in  the  attempts  to  grow  single  crystals 
of  the  complexes  supplied,  spectra  were  taken  of  the  complexes  in 
their  powder  or  microcrystalline  form.   The  resulting  spectra  all 
exhibited  extremely  broad  lines  at  room  temperature.   Initially  attri- 
buting the  broad  lines  simply  to  spin- lattice  relaxation  effects,  a 
series  of  experiments  was  performed  on  four  particular  complexes  which 
exhibited  anisotropic  line  shapes  at  room  temperature.   These  experi- 
ments revealed  temperature  dependent  behavior  of  the  spectra  of  all 
four  complexes.   Figure  4  shows  the  variation  of  the  powder  spectra  of 
the  four  different  tetrazole  complexes  with  temperature.   The  four  com- 
plexes studied  were: 

Figure  4A: 

Figure  4B 

Figure  4C :   1,5-pentamethy lene  tetrazole  with  CuBr 

Figure  4D:   l-benzyl-5-amino  tetrazole  with  Cu  Br 

Encouraged  by  the  fact  that  for  some  of  the  complexes,  the  spectral 
lines  narrowed  at  lower  temperature,  it  was  decided  to  try  work  at 
liquid  helium  temperature.   Initial  preparations  with  a  spectrometer 
similar  to  that  used  for  the  liquid  nitrogen  work  proceeded  slowly  and 


l-methyl-5-ethyl  tetrazole  with  CuCl 
l-methyl-5-amino  tetrazole  with  CuCl 
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FIGURE:     4A 
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FIGURE       4  B 
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FIGURE  4D 
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with  much  difficulty,  but  finally  a  powder  spectrum  of  l-methyl-5-amino 
tetrazole  copper  (II)  chloride  complex  was  obtained  at  liquid  helium 
temperature.   The  result  is  shown  in  Fig.  4A.   It  was  concluded  on  the 
basis  of  this  experiment  that  the  liquid  helium  work  would  bear  little 
fruitful  results,  and  so  efforts  in  this  direction  were  terminated. 

Room  temperature  powder  spectra  of  the  complexes  which  exhibited 
isotropic  line  shapes  were  taken  and  the  line  width  measured.   Room 
temperature  powder  spectra  of  the  three  copper  (I)  complexes  were  also 
taken.   The  results  are  listed  in  Table  2. 

Although  no  positive  explanation  can  be  given  for  the  temperature 
dependent  behavior  of  some  of  the  complexes,  hypotheses  may  be  advanced 
to  explain  the  effects.   Typical  behavior  is  exemplified  by  l-methyl-5- 
ethyl  tetrazole  with  CuCl   in  Fig.  4A.   The  spectrum  seems  to  alter- 
nately broaden  and  narrow  as  the  temperature  is  raised  or  lowered.   The 
importance  of  some  line  broadening  effects  will  now  be  considered  to 
explain  this  behavior. 

At  high  temperature,  the  effects  of  spin  lattice  relaxation  are 
important.   In  addition,  there  exists  the  possibility  of  rapid  exchange 
of  the  tetrazole  ligand  between  two  positions  it  could  occupy  (Fig.  1, 
top)  in  the  puckered  polymeric  chain.   This  possibility  is  advanced  in 
light  of  the  discovery  of  previous  workers  [9]  that  the  tetrazole  ligand 
is  in  general  a  poor  donor,  and  therefore  does  not  interact  strongly 
with  the  central  copper.   The  polymeric  chain  would  then  be  free  to 
dynamically  "repucker"  by  making  and  breaking  bonds  with  the  tetrazole. 
Spin- lattice  effects  would  tend  to  broaden  the  spectral  line  and  ex- 
change tend  to  narrow  it.   At  any  particular  temperature  the  dominance 
of  one  effect  or  the  other  would  depend  upon  the  structure  of  the 
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complex  and  the  temperature  dependence  of  the  exchange  mechanism.   At 
lower  temperatures,  where  the  spin-lattice  effects  are  greatly  diminished 
and  the  lattice  vibrations  are  very  small,  it  is  possible  for  the  copper 
atoms  to  be  close  enough  so  that  the  electron  is  affected  by  both  copper 
nuclei.   This  effect  leads  to  dipolar  relaxation  and  line  broadening. 

C.   SINGLE  CRYSTAL  WORK 

Of  the  many  copper  halide  tetrazole  complexes  supplied  by  R.  A. 
Henry  of  China  Lake,  only  the  2-methyl-5-amino  copper  (II)  chloride  and 
bromide  tetrazole  complexes  were  in  the  form  of  crystalline  chips  large 
enough  to  be  treated  as  single  crystals.  Numerous   attempts  at  growing 
single  crystals  of  the  other  complexes  were  made,  and  all  but  one  at- 
tempt met  with  failure.   The  conventional  method  of  growing  crystals  by 
slow  evaporation  from  solution  was  impractical  because  the  copper 
halide  tetrazole  complexes  which  were  supplied  were  insoluble  in  all 
common  organic  solvents  tried.   Basically,  all  of  the  techniques  for 
growing  single  crystals  were  based  on  the  method  of  slow  diffusion  of 
the  copper  halide  and  the  tetrazole  dissolved  in  a  suitable  solvent. 
The  slow  diffusion  technique  was  necessary  because  of  the  tendency  of 
the  complexes  to  immediately  precipitate  when  a  solution  of  tetrazole  in 
solvent  was  poured  into  a  solution  of  copper  halide  in  the  same  solvent. 
The  complexes  have  a  very  large  formation  constant  [7] . 

The  successful  attempt  at  growing  a  single  crystal  of  2-methyl-5- 
amino  copper  chloride  tetrazole  complex  yielded  a  crystal  approximately 
ten  times  the  size  of  any  of  the  small  crystals  initially  supplied  from 
China  Lake.   A  few  grams  of  copper  chloride  was  placed  in  the  bottom  of 
a  long  pyrex  tube  sealed  at  one  end.   The  tube  was  filled  with  glass 
wool  and  a  few  grams  of  the  2-methy 1-5-amino  tetrazole  placed  on  top  of 
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the  column  of  glass  wool.   The  entire  tube  was  filled  with  absolute 
ethanol  and  then  stoppered  and  allowed  to  stand  at  room  temperature. 
After  a  period  of  two  months,  a  crystal  about  the  size  of  a  pea  was  ob- 
served to  be  clinging  to  the  column  of  glass  wool.   Due  to  the  time 
period  involved  in  this  technique,  it  was  not  tried  with  other  tetra- 
zoles.   Of  all  the  techniques  for  growing  the  single  crystals  of  copper 
halide  tetrazoles,  this  one  showed  the  most  promise.   However,  the  com- 
plex from  which  the  crystal  was  grown  had  already  grown  smaller  single 
crystals.   This  was  not  the  case  for  the  rest  of  the  complexes  which 
came  in  powder  or  microcrystalline  form. 

The  results  of  the  single  crystal  work  with  the  two  complexes  men- 
tioned above  is  reported  in  Table  3.   The  spectra  of  the  single  crystals 
showed  one  sharp,  isotropic  line  at  any  orientation  of  the  magnetic 
field  which  means  that  there  exists  only  one  magnetically  distinct  site 
in  the  crystal.   In  both  crystals,  the  g  values  are  all  different  with 
two  out  of  each  set  differing  only  by  0.03.   Both  of  these  observations 
are  consistent  with  the  proposed  polymeric  structure  of  a  2-subs tituted 
tetrazole  (Fig.  1,  bottom).   The  two  g  values  which  are  almost  equal 
correspond  to  the  g  values  taken  when  the  magnetic  field  is  parallel 
to  or  perpendicular  with  the  plane  of  the  polymer.   The  proposed  struct- 
ure also  has  only  one  magnetically  distinct  site. 

The  g  values  may  be  used  to  make  a  rough  calculation  as  to  the 
amount  of  covalency  in  the  copper- tetrazole  bond.   Assuming  the  struct- 
ure at  the  bottom  of  Fig.  1  to  be  the  structure  of  the  complexes  under 
consideration  here,  two  formulae  for  the  g  values  of  an  octahedral 
copper  complex  with  a  tetragonal  distortion  corresponding  to  elonga- 
tion along  a  four- fold  symmetry  axis  may  be  developed.   The  development 
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is  similar  to  the  development  of  the  g  value  formulae  derived  in  the 
introduction.   The  final  formulae  are: 


g   =  2.002  +|£  P 
g   =  2.002  +  7*  P 


5  =  the  spin  orbit  constant  for 
copper  =  829  cm 

P  =  the  one  electron  popula- 
tion of  the  copper  d 
orbital  [2] 


Jonassen  [5]   reported  the  values  for  A   for  both  complexes  but  gave  no 
value  for  A  .  With  the  above  formulae,  P  may  be  calculated  and  then 
used  to  predict  A  .   The  results  are  given  in  Table  3. 
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Position  and  kind 

of  substituent  on      Composition        g  values 
tetrazole 


linewidths 
(gauss) 


2-methyl-5  amino     2(C  H  N  )  CuCl     2.043  +  0.002    24.5  +  0.5 

2.073  29.4 


2.266  49.0 


2-methyl-5   amino  2(C   H  N,.)    CuBr  2.043  +   0.002  22.0  +  0.5 

2.073  34.3 

2.218  36.8 


(g        -    2.002)A 

Now   use  P   =  — — (1) 

8^ 


In   the   computation  of  A.    = *'_    0„„N  (2) 

2    (g   -  2.002) 


2%? 
'1 


Chloride  Complex  Bromide  Complex 

Given:   g    =  2.266  g^   =  2.218 

A   =  14100  cm"1  A   =  14380  cm"1 

Obtain  from  (1)  P  =  0.56  P   =  0.47 


Using  these  values  of  P  and  taking  g,  =  average  of  g   and  g 

°   1  xx      yy 

=  2.058 
Result:   A   =  16580  cm"1  A   =  13920  cm" 


TABLE  3.   RESULTS  OF  SINGLE  CRYSTAL  WORK 
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VII.   THEORY  OF  THE  ESR  OF  TRAPPED  ORGANIC  RADICALS 

The  effective  spin  Hamiltonian  describing  the  interaction  of  the 
nuclei  and  electron  of  a  trapped  organic  radical  with  the  magnetic 
field  and  with  one  another  has  previously  been  written  down  in  the 
section  on  the  ESR  of  transition  metals.   It  will  be  listed  again  to 

aid  in  the  following  discussion: 

— ♦  — > 

K  =  $H  •  g  •  S-Lv,M  -I.  +SS  -A  •  1 

i   i      i    i  i 

In  the  section  on  transition  metals,  the  first  term  of  the  effective 
spin  Hamiltonian  was  thoroughly  described.   For  trapped  organic  radi- 
cals, the  g  tensor  of  this  term  deviates  only  slightly  from  the  free 
spin  g  value.   It  is  symmetrical  and  its  principal  values  are  always 
close  to  the  free  spin  g  value.   The  second  term  represents  the  inter- 
action of  the  nuclear  spins  with  the  magnetic  field.   The  effect  of 
this  term  is  primarily  to  relax  the  selection  rule  Am  =  0. 

Reference  10  contains  an  excellent  description  of  the  third  term 
which  represents  the  hyperfine  interaction  between  the  electron  spin  S 
and  the  nuclear  spin  I  described  by  the  second-rank  tensor  A.  Insofar 
as  g  is  nearly  isotropic,  the  hyperfine  tensor  A  will  be  a  symmetrical 
3X3  matrix  which  can  be  diagonalized  and  resolved  into  isotropic  and 
anisotropic  components.  The  effective  Hamiltonian  H1  for  the  hyper- 
fine interaction  energy  can  be  written: 


K.  .  ftrgH   t2(Q)  -  .  -  .  gEW 


— »     — »   _  Q        — t     — •   — »     — ♦    -S 

S  •  Ir   -  3(S  •  r)(I  •  r)r 


The  eigenvalues  of  H'  are  the  principal  values  of  A,  and  the  isotropic 
and  anisotropic  components  of  A  correlate  with  the  first  and  second 
terms,  respectively. 
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The  first  term  in  the  equation  for  H'  gives  rise  to  a  nonzero  hyper- 
fine  splitting  only  if  there  is  a  finite  probability  of  finding  the  un- 
paired electron  at  the  interacting  nucleus,  as  indicated  by  the  factor 

2 
ill  (0).   For  this  reason  the  isotropic  interaction  is  usually  called  the 

contact  or  Fermi  interaction.  It  is  apparent  that  the  isotropic  hyper- 
fine  interaction  A  will  only  be  nonzero  if  ilr(o)  is  nonzero;  that  is,  if 
the  orbital  of  the  unpaired  electron  has  s-character. 

The  second  term  in  the  equation  for  H1  imparts  to  the  effective 
Hamiltonian  an  anisotropic  or  direction  dependent  character.   This 
anisotropic  part  of  the  hyperfine  tensor  is  a  traceless  tensor  which 
describes  the  hyperfine  interaction  energy  consequent  upon  the  classi- 
cal interaction  of  the  electronic  and  nuclear  dipoles. 

In  general,  the  hyperfine  interaction  of  n  nuclei  of  spin  1/2  will 
give  rise  to  a  spectrum  consisting  of  2   equally  intense  lines.   If  two 
or  more  equivalent  spin  1/2  nuclei  interact,  superposition  of  some  of 
the  lines  results  in  a  spectrum  of  n  +  1  lines  of  binomial  intensity 
distribution. 
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VIII.  BACKGROUND  ON  THE  ESR  SPECTRA  OF  TRAPPED  ORGANIC  RADICALS 

A  considerable  amount  of  work  has  been  done  on  the  subject  of  the 
ESR  spectra  of  radiation  damaged  organic  crystals.   Naturally,  the  work 
has  covered  those  materials  from  which  single  crystals  are  easily  ob- 
tained.  Crystals  of  saturated  carboxylic  acids  and  their  salts,  the 
simpler  amino  acids  and  certain  unsaturated  acids,  when  irradiated  with 
x-rays,  gamma  rays,  or  high  energy  electrons,  have  been  found  to  pro- 
vide suitable  hosts  for  the  trapped  radicals  produced  by  such  irradia- 
tion.  Analysis  of  the  spectra  from  radiation  damaged  crystals  of  these 
materials  is  primarily  concerned  with  the  identification  of  the  trapped 
radical(s)  via  the  interpretation  of  the  proton  hyperfine  interaction. 
One  type  of  radical  which  is  commonly  detected  is  a  rr  -  electron  radi- 
cal centered  on  a  carbon  atom.   In  such  a  radical  the  unpaired  elect- 
ron predominantly  occupies  a  carbon  2p  orbital  perpendicular  to  the 

i 
plane  of  the  radical. 

Kispert  and  Rogers  [13]  reported  the  interpretation  of  the  ESR 
spectrum  of  single  crystals  of  sodium  acetate  trihydrate  irradiated 
with  1  Mev  electrons  at  77   K.   At  this  temperature,  they  observed  a 
four  line  pattern  with  intensity  ratio  of  1:3:3:1  and  proton  hyper- 
fine splitting  of  22.5  gauss.   Each  of  these  lines  had  two  satellites 
spaced  5.6  gauss  on  either  side  to  form  a  twelve  line  pattern.   They 
interpreted  the  main  four  line  pattern  to  be  due  to  the  methyl  radical, 
•CH  .   The  satellite  lines  were  attributed  to  spin  flip  transitions 
arising  from  a  simultaneous  change  in  the  spin  state  of  the  odd  elect- 
ron on  the  methyl  radical  and  of  the  neighboring  protons  of  the 
molecules  of  water  of  crystallization.   One  of  the  first  workers  to 
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investigate  the  ESR  of  an  x-irradiated  single  crystal  of  methyl  malonic 
acid  was  Heller  [14].   He  irradiated  the  malonic  acid  at  room  tempera- 
ture and  observed  the  spectra  at  room  temperature,  77  K,  and  4  K. 
Heller  interpreted  his  spectra  in  terms  of  the  two  species  CH_C(C00H) 
and  CH„CHC00H.   He  concentrated  his  efforts  on  the  first  and  con- 
cluded, from  the  spectra,  that  the  three  methyl  protons  are  equivalent, 
and  the  principal  values  of  their  hyperfine  coupling  tensor  are  75.4 
(along  the  C-CH»  bond),  68.8  and  68.6.   He  also  concluded  that  the 
isotropic  component  of  the  coupling  tensor  (and,  therefore,  the  spin 
density  on  the  methyl  protons)  is  positive,  and  that  the  methyl  group 
executes  nearly  free  rotation  about  the  C-C  bond  at  4.2  K.   Tamura, 
Collins,  and  Whiffen  [15]  irradiated  various  malonic  acids  at  77   K 
and  observed  the  ESR  spectra  as  the  temperature  was  raised.   They 
found  that  the  radical  formed  at  77°  K  was  of  the  form  R-CH  -CH(COOH) , 
and  that  as  the  temperature  was  raised,  a  reaction  took  place  in  which 
the  original  radical  disappeared  and  a  radical  of  the  form  RCH  C(COOH) 
showed  up  in  its  place. 

J.  R.  Morton  [161  reported  the  ESR  spectrum  of  (CH  )  CO  "   trapped 
in  a-aminoisobutyric  acid  at  40   K.   He  concluded  that  at  this  tempera- 
ture, one  of  the  methyl  groups  is  undergoing  virtually  free  rotation 
about  the  C-CH   bond,  while  the  other  methyl  group  is  undergoing  hinder- 
ed rotation.   Box  and  Freund  [17]  investigated  x-irradiated  a-amino- 
isobutyric  acid  irradiated  at  77  K.   They  allowed  the  crystal  to  warm 
gradually  and  observed  the  resulting  spectra.   They  detected  distinct 
intermediate  ESR  spectra  due  to  preliminary  conformations  of  the  radi- 
cal present  at  room  temperature.   They  found  the  intermediate  spectra 
to  be  better  resolved  and  easier  to  interpret  than  the.  final  spectrum. 
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They  also  found  that  cooling  the  intermediate  fragment  to  77   K  showed 
no  change  in  the  spectrum,  but  that  cooling  the  room  temperature  frag- 
ment did  show  a  changing  spectrum.   They  interpreted  this  behavior  as 
resulting  from  different  conformations  of  (CH  )  CCOOH. 

Miyagawa  and  Gordy  [18]  first  reported  the  results  of  their  investi- 
gation of  room  temperature  irradiated  single  crystals  of  alanine.   They 
proved  that  the  stable  free  radical  produced  by  the  irradiation  was  of 
the  form  CH„CHR,  where  R  is  a  group  which  has  no  nuclei  with  detectable 
coupling.   The  hydrogens  of  the  methyl  group  of  the  radical  were  shown 
to  have  equivalent,  isotropic  coupling  of  26  gauss  each,  essentially 
independent  of  the  frequency  of  observation.   Horsfield,  Morton,  and 
Whiffen  [19  and  20]  observed  the  change  in  the  ESR  spectrum  of  the  room 
temperature  gamma  irradiated  single  crystal  of  1-a-alanine  on  cooling 
to  77   K,  and  interpreted  their  results  in  terms  of  the  cessation  of 
rotation  of  the  methyl  group  of  the  trapped  radical  CH  CHC00H.   The 
isotropic  hyperfine  couplings  of  the  three  hydrogens  of  the  methyl 
group  were  found  to  be  120,  76,  and  14  Mc/sec  compared  to  the  value  of 
70  Mc/sec  each  at  300  K.   In  the  intermediate  temperature  range  bet- 
ween 100   K  and  200   K,  the  spectra  were  found  to  differ  from  the 
spectra  above  and  below  these  temperatures.   From  the  line  width  change 
near  100  K  and  200  K,  an  estimate  of  approximately  4  kcals/mole  was 
made  for  the  activation  energy  hindering  the  rotation  of  the  methyl 
group.   Miyagawa  and  Itoh  [21]  determined  that  room  temperature  gamma 
irradiated  single  crystlas  of  dl-alanine  grown  from  heavy  water  ex- 
hibited the  same  hyperfine  structure  as  those  crystals  grown  from  nor- 
mal water  except  that  the  line  widths  were  considerably  sharper.   They 
observed  the  spectra  as  the  temperature  was  lowered  and  concluded  that 
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the  rotation  of  the  methyl  group  was  hindered  at  lower  temperatures. 
They  used  a  modified  Bloch  equation  to  estimate  the  frequency  of  rota- 
tion of  the  methyl  group  from  the  absorption  patterns  at  various  temper- 
atures, and  by  applying  the  classical  theory  of  rate  processes  to 
analyse  the  relation  between  the  frequency  of  rotation  and  the  tempera- 
ture, they  found  the  potential  barrier  to  internal  rotation  of  the  methyl 
group  to  be  3.6  kcal/mole. 

Sinclair  and  Hanna  [22  and  23]  examined  the  ESR  spectra  of  1-alanine- 

13  o 

1-   C,  x-irradiated  at  80  K.   They  concluded  that  the  low  temperature 

initial  species  has  a  high  unpaired  electron  population  on  the  carboxyl 

carbon  and  that  as  the  low  temperature  radical  is  allowed  to  decay 

thermally  to  the  fragment  CH  CHCOOH  at  140-150°  K,  there  is  a  transfer 

of  the  unpaired  electron  from  the  carboxyl  group  to  the  a-carbon.   The 

CH~CHC00H  is  at  a  different  orientation  in  the  crystal  than  the 

CH  CHCOOH  formed  by  irradiation  at  room  temperature.   Upon  warming  to 

room  temperature,  a  reorientation  of  the  fragment  occured  which  gave 

the  usual  room  temperature  ESR  spectrum.   A  feature  of  the  intermediate 

fragment  was  that  the  methyl  protons  remained  equivalent  even  when  re- 

cooled  to  80  K.   This  was  in  contrast  to  the  room  temperature  fragment 

where  the  rapid  interconversion  of  methyl  protons  was  frozen  out  by 

o 
cooling  to  77   K. 

Wells  and  Box  [24]  studied  the  rotation  of  the  methyl  group  in  the 

free  radical  CjfcCH  -C0NH-C(CH  )C00H  by  ESR  and  ENDOR  spectroscopy.   They 

analyzed  their  results  using  a  quantum  mechanical  hindered  rotor  model 

with  Mathieu  functions  as  basis  functions.   Their  experimental  results 

indicated  a  barrier  height  of  less  than  0.5  kcal/mole  for  the  methyl 

group  rotation. 
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In  recent  paper  on  x- irradiated  single  crystals  of  zinc  acetate 

dihydrate,  Tolles,  Crawford,  and  Valenti  [25]  reported  the  presence 

of  the  methyl  radical  in  the  case  of  irradiation  at  77   K,  with  the 

subsequent  disappearance  of  this  radical  and  the  appearance  of  *CH  CO 

as  the  temperature  was  raised.   They  proposed  a  mechanism  to  explain 

the  decay  of  the  first  species  and  corresponding  increase  in  the 

second  species.   The  mechanism  involves  abstraction  of  a  hydrogen  atom 

by  the  methyl  radical  from  a  neighboring  acetate  ion.   They  supported 

this  hypothesis  by  synthesizing  Zn(CH„DCO.)  • 2H„0     ,   ,     ,   .     , 
Jr  2   2^2   and  showed  that  the 

decay  and  growth  of  the  above  species  exhibits  an  isotope  effect. 

Spectra  taken  over  a  range  of  temperatures  revealed  both  the  slow  and 

rapid  exchange  of  the  hydrogen  atoms  seen  previously  [e.g.  16,17,18], 

with  the  barrier  to  internal  rotation  of  the  -CH_  group  being  5  kcal/ 

mole. 
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IX.      EXPERIMENTAL  WORK 

A.  GROWING  OF  CRYSTALS 

Single  crystals  of  strontium  acetate  hemihydrate  were  grown  by  slow 
evaporation  from  a  saturated  aqueous  solution  at  room  temperature. 
Crystals  in  which  heavy  water  was  exchanged  for  the  water  of  hydration 
were  grown  by  making  a  saturated  solution  of  the  strontium  acetate 
hemihydrate  in  997o  DO  at  room  temperature,  and  then  allowing  slow 
evaporation  of  the  solution.   The  single  crystals  crystallized  in  the 
form  shown  in  Fig.  5  which  enabled  definition  of  an  orthogonal  crystal 
axis  system.   The  crystals  cleave  easily  parallel  to  the  XY  plane. 
The  XYZ  crystal  axis  system  is  shown  in  Fig.  5.   The  numbers  along  the 
edges  of  the  crystal  are  the  angles  the  crystal  faces  make  with  the 
XY  plane. 

B.  IRRADIATION  OF  THE  CRYSTAL 

The  single  crystals  were  exposed  to  x-irradiation  at  50  KV  and  40 
ma  from  a  copper  target.   The  crystals  irradiated  at  room  temperature 
were  placed  on  scotch  tape  and  mounted  in  the  open  windows  of  the  x-ray 
machine  allowing  the  port  covers  to  be  closed  at  all  times.   For  room 
temperature  irradiations,  the  crystals  were  exposed  for  one  hour.   For 
irradiation  at  liquid  nitrogen  temperature,  the  crystals  were  placed  in 
a  quartz  tube  sealed  at  one  end.   The  other  end  was  stoppered  to  pre- 
vent condensation  of  liquid  air  into  the  tube.   The  tube  with  crystals 
was  placed  in  a  styrofoam  ice  bucket  and  precooled  gradually  by  slowly 
pouring  liquid  nitrogen  into  the  bucket  and  allowing  the  cold  gas  to 
flow  over  the  tube.   Precooling  was  necessary  to  prevent  shattering  of 
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FIGURE  5.   A  single  crystal  of  Sr(CH_CO_)   ■  %H.O  showing  the  angles 

5       2-2  2 

between  the  edges  and  faces  of  the  crystal,  and  the  ortho- 
gonal crystal  axis  system. 
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the  crystal  when  exposed  to  a  sudden  temperature  change.   Eventually, 
the  tube  was  immersed  completely  in  liquid  nitrogen.   The  bucket  was 
mounted  on  a  ring  stand  with  a  lead  shield  attached,  and  a  wooden  cover 
with  a  hole  drilled  in  the  top  of  the  bucket.   The  bucket  was  placed 
in  front  of  the  x-ray  port,  and  the  tube  positioned  in  the  x-ray  beam 
so  that  the  crystal  received  maximum  dosage  for  a  two  hour  period.   The 
position  of  the  tube  in  the  beam  was  checked  by  darkening  the  room  and 
lowering  a  piece  of  fluorescent  paper  into  the  bucket.   The  tube  and 
crystal  formed  a  shadow  on  the  fluorescent  paper. 

C.   CRYSTAL  ALIGNMENT 

The  orientation  of  the  crystal  in  the  spectrometer  cavity  had  to 
be  determined  as  precisely  as  possible  in  order  to  determine  the  loca- 
tion of  the  magnetic  field  with  respect  to  the  crystal  axis  system. 
The  crystal  was  mounted  on  a  quartz  rod  which  had  two  flats  ground 
parallel  to  each  other  on  both  ends  of  the  rod.   The  crystal  was  held 
on  the  rod  by  means  of  a  platinum  spring  attached  to  the  rod  by  a 
piece  of  scotch  tape.   The  end  opposite  the  mounted  crystal  was  clip- 
ped to  an  angle  indicator  positioned  over  the  center  of  the  spectro- 
meter cavity.   The  crystal  was  positioned  on  the  rod  by  use  of  a  micro- 
scope, one  of  the  eyepieces  of  which  had  crosshairs.   The  orientations 
of  the  crystal  on  the  rod  were  chosen  so  that  the  well  defined  edges  of 
crystal  were  parallel  or  perpendicular  to  the  edges  of  the  rod.   The 
angles  between  the  faces  and  edges  of  the  crystal  were  determined  with 
an  optical  goniometer.   For  liquid  nitrogen  work,  the  crystal  alignment 
was  accomplished  without  warmup  of  the  crystal  by  keeping  the  crystal 
immersed  in  a  bath  of  liquid  nitrogen  under  the  microscope. 
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D.   TAKING  THE  SPECTRA 

The  spectra  were  taken  with  the  spectrometer  previously  described 
for  the  liquid  nitrogen  work  on  copper  tetrazole.   The  temperature  was 
controlled  with  the  temperature  control  unit  attachment  to  the  spectro- 
meter used  for  the  copper  tetrazole  work. 
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X.   RESULTS 

A.   ROOM  TEMPERATURE  SPECIES 

Spectra  obtained  from  x-ray  damaged  single  crystals  of  strontium 
acetate  hemihydrate  at  room  temperature  were  found  to  depend  upon  the 
history  of  the  crystal  as  well  as  the  orientation  of  the  magnetic  field 
relative  to  the  crystal  axis  system.  As  can  be  seen  in  Fig.  6  (top), 
crystals  irradiated  at  liquid  nitrogen  temperature,  slowly  warmed  to 
room  temperature,  and  then  observed  at  room  temperature  did  not  exhibit 
the  extra  central  line  in  the  spectrum  as  did  the  crystals  irradiated 
and  observed  at  room  temperature  (Fig.  6,  bottom).   Both  spectra  were 
obtained  at  the  same  crystal  orientation  of  &   =  0  ,  9  =  90  .   Reorient- 
ation of  the  magnetic  field  resulted  in  a  coalescence  of  the  two 
1:3:3:1  patterns  into  what  appears  to  be  a  1:4:6:4:1  pattern  as  in 
Fig.  7.   The  orientations  for  the  magnetic  field  are  from  top  to  bottom, 
respectively  4  =   0°,  9  =  45° ;  i  =   0°,  9  =  80° ;  4   =   0°,  9  =  90°.   The 
spectra  can  be  interpreted  in  terms  of  a  free  radical  species  with 
three  equivalent  hydrogens  split  by  a  fourth,  inequivalent  hydrogen. 
The  spectra  obtained  also  indicate  that  there  is  only  one  magnetically 
distinct  site  in  the  crystal  for  the  room  temperature  species. 

For  some  orientations  of  the  magnetic  field,  the  base  line  of  the 
spectrum  was  not  straight,  but  rather  convex  upward.   The  sweep  range 
was  therefore  changed  from  250  gauss  to  1000  gauss  and  the  spectra  in 
Fig.  8  were  obtained.   The  spectra  indicate  the  possible  presence  of 
a  triplet. 
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FIGURE     6 
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figure:     t 
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e  =  60 


100  G 


%  =   ou 

e  =  120 


FIGURE  8. 


Broad  line  species  seen  in  SrCCH^COp 
observed  at  room  temperature. 
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\   HO  irradiated  and 


Several  spectra  of  the  room  temperature  species  were  taken  at 
different  magnetic  field  orientations.   The  coupling  constant  for  the 
inequivalent  hydrogen  was  not  always  obtainable  due  to  the  overlapping 
of  the  spectral  lines  at  several  orientations.   However,  the  coupling 
constant  for  the  CH„  group  was  always  obtainable.   The  methyl  group 
coupling  constants  and  corresponding  magnetic  field  orientations  were 
plotted  on  a  s tereographic  net,  and  an  estimate  of  the  positions  of 
the  principal  axis  system  and  the  coupling  constants  along  the  princi- 
pal axis  was  made,  assuming  axial  symmetry  for  the  rotating  methyl 
group.   The  results  are: 

A   =  27.1  +  0.4  gauss,  direction  cosines  (-0.32,  0.38,-0.87) 

A   =  24.2  +  0.4  gauss 

1 

The  errors  in  the  direction  of  the  principal  axis  are  estimated  to  be 
ten  degrees  and  for  the  coupling  constants,  0.5  gauss.   The  coupling 
constants  and  magnetic  field  orientations  for  the  inequivalent  hydrogen 
were  least  squared  [26].   The  results  are,  for  the  three  principal 
values  of  the  A  tensor: 

(a)  11.61  +  0.40  gauss 

(b)  32.98 

(c)  19.32 

The  direction  cosines  for  the  principal  axes  are: 

(a)  (0.72,  0.69,  0.80) 

(b)  (-0.18,  0.73,  0.98) 

(c)  (0.67,  -0.72,  0.18) 

in  the  XYZ  axis  system  in  Fig.  5.   The  data  used  is  in  Table  4.   Figure 
9  shows  the  angular  dependence  of  the  coupling  constant.   The  graphs 
indicate  that  the  methyl  group  coupling  constants  show  a  cos9  dependence, 
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9  a   3  (gauss) 


0  25  25.54  +  0.35 

30  25.24 

35  25.24 

40  25.06 

45  24.90 

50  24.56 

55  24.56 

60  24.38 

65  24.38 

70  24.20 

80  24.38 


70  30  25.58 

45  24.90 

60  24.38 

75  24.38 

90  24.20 

105  24.90 

120  25.41 

135  26.28 

150  26.45 


45  0  26.10 

15  25.93 

30  25.54 

45  25.24 


TABLE  4.   Coupling  constants  and  magnetic  field  orientations  for  the 
room  temperature  species. 
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Table    4   continued, 


(45) 


e 


60 
75 
90 
105 
120 
135 
150 


CHT  (  s 

a   3  (gauss) 

24.55  +  0.35 

24.38 

24.20 

24.20 

25.24 

25.76 

26.28 


90 


0 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 


26.28 
25.93 
25.24 
25.06 
24.90 
24.38 
24.55 
25.06 
25.59 
26.10 
26.45 
26.23 


25 


30 

45 

60 

75 

90 

105 

120 

135 

150 

165 


25.42 
24.72 
24.38 
24.38 
24.20 
25.06 
25.42 
26.28 
26.45 
26.80 
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Table  4  continued. 

_J_  9 


a      3    (gauss) 

25, 

.93  + 

0, 

,35 

25, 

.93 

25, 

,54 

25, 

,06 

24, 

.38 

24, 

,77 

24, 

.56 

24, 

,90 

25, 

.93 

26, 

.45 

26, 

.97 

26, 

.62 

115  0 

15 
30 

45 

60 

75 

90 

105 

120 

135 

150 

165 


<j>  9  a    (gauss) 

0  25  17.12  +  0.35 

30  15.04 

35  13.83 

40  12.10 

45  11.24 

50  11.06 

55  11.93 

60  13.83 

65  15.73 

70  16.94 

25  30  15.22 

45  10.54 

60  12.44 

75  17.29 

120  29.39 

135  30.60 

150  30.26 
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Table   4   continued, 


<t>  8  a        (gauss) 

45  15  20.74  +  0.35 

30  16.94 

45  13.31 

60  12.80 

75  15.90 

90  20.22 

70  30  20.22 

45  18.15 

60  16.94 

75  17.46 

90  14.01 
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a   m  gauss 
H    6 


9  in  degrees 


FIGURE  9A. 


H  atom  coupling  constant  versus  orientation  of  magnetic 
field  in  XYZ  axis  system.  0=0° 
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FIGURE  9B.  <f>   =   25 
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a   in  gauss 
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FIGURE  9C.  p   =  45 


6  in  degree  s 
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FIGURE  9D.  <t>   =  70 


8  in  degrees 
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aCH    in  gauss 
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FIGURE  9E.   CH_  group  coupling  constant  versus  orientation  of 


magnetic  field  in  XYZ  axis  system.  <j>   =   0 


=  no 
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a     in  gauss 
CH3 
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FIGURE  9F.  <t>   =  25 
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a     in  gauss 
CH3 
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FIGURE  9G.  i   *  45 
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FIGURE  9H.  <t>   =  7  0 
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aCH   in  gauss 
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FIGURE  91.  i   »  90 
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HCH    in  gauss 
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FIGURE  9J.  <h   =  L15 
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Powdered  strontium  acetate  hemihydrate  was  recrystallized  from  997o 
DO  and  the  resulting  radiation  damaged  crystals  showed  the  same  spect- 
rum as  single  crystals  recrystallized  from  normal  water.   It  may  be  con- 
cluded that  the  water  of  hydration  does  not  enter  into  reactions  involv- 
ing the  initial  free  radical  fragments. 

An  interesting  sidelight  to  the  room  temperature  work  is  that  a 
crystal  of  zinc  acetate  dihydrate  which  had  been  x-ray  damaged  and  then 
set  aside  for  one  month  was  looked  at  on  the  spectrometer.   Surpris- 
ingly, the  spectrum  obtained  was  the  same  1:4:6:4:1  pattern  obtained 
for  strontium  acetate  hemihydrate.   Apparently,  the  radical  species 
identified  by  Tolles,  Crawford,  and  Valenti  [25]  was  an  intermediate 
in  a  reaction  with  a  very  small  rate  constant  or  large  activation 
energy  which  leads  to  a  final  product  similar  to  that  in  strontium 
acetate.   The  spectrum  of  the  species  in  zinc  acetate  dihydrate  is 
shown  in  Fig.  10. 

B.   INTERMEDIATE  SPECIES 

Single  crystals  of  strontium  acetate  hemihydrate  irradiated  at 
liquid  nitrogen  temperature  and  observed  at  -100  C  exhibited  a  1:3:3:1 
pattern  which  split  into  eight  lines  of  equal  intensity  as  the  tempera- 
ture was  lowered  to  -160  C.   Spectra  for  the  species  present  at  -100  C 
were  taken  over  several  orientations  of  the  magnetic  field,  the  coupling 
constants  and  magnetic  field  orientations  plotted  on  a  stereographic 
net,  and,  assuming  axial  symmetry  for  a  rotating  methyl  group,  the 
hyperfine  coupling  constants  and  direction  cosines  of  the  species  rela- 
tive to  the  XYZ  axis  system  in  Fig.  5  were  calculated  to  be  [26]: 
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A  =  12.84  +  0.30  gauss 

A  =  10.04  +  0.30  gauss 

X  =  (0.81,  0.58,  0.06) 

Y  =  (-0.59,  0.80,  0.12) 

Z  =  (0.02,  -0.14,  0.99) 

The  data  fed  into  the  program  is  in  Table  5.   The  average  isotropic 
hyperfine  coupling  constant  is  9.2  gauss. 
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JL_  6  aCH3(Mc/sec) 

0  0  35.89+1.00 

20  36.83 

40  36.27 

60  32.16 

70  31.78 

80  30.28 

90  29.35 

100  27.86 

110  28.60 

120  28.42 

140  31.03 

160  33.84 

90  0  35.33 

20  36.36 

40  34.96 

60  31.03 

70  29.16 

80  28.42 

90  27.86 

100  27.29 

110  28.04 

120  28.98 

140  31.22 

160  34.21 


TABLE  5.   Coupling  constants  and  magnetic  field  orientations  for  the 
-100°  C  species. 
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Table  5  continued 
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0 
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CH 
a   3  (Mc/sec) 


35 

.71  +  0.35 

34 

.40 

31 

.41 

28 

.79 

27 

.48 

27 

.11 

28 

.60 

29 

.72 

31 

.59 

33 

.46 

35 

.33 

35 

.71 

36 

.64 

35 

.33 

33 

28 

30 

.10 

29 

54 

28 

23 

28. 

04 

27. 

86 

28. 

60 

28. 

60 

32. 

90 
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The  spectra  of  this  same  species  were  obtained  as  a  function  of 
temperature.   The  computer  program  following  the  appendix,  which  is 
based  upon  the  mathematical  model  for  motional  and  exchange  narrowing 
of  spectral  lines,  was  used  to  predict  the  spectrum  of  the  rotating 
methyl  group  for  various  rotational  frequencies.   In  the  drawing  be- 
low, the  eight  line  pattern  resulting  from  splitting  by  three  equiva- 
lent protons  corresponds  to  the  zero  rotation  case.   In  the  case  of 
rotation  of  the  methyl  group,  all  of  the  lines  except  the  two  outside 
lines  exchange,  and  the  extent  to  which  they  exchange  is  a  function  of 
the  frequency  at  which  the  methyl  group  is  rotating: 


aaot 


Various  rotational  frequencies  were  fed  into  the  above  mentioned  com- 
puter program,  and  the  computed  spectra  compared  with  the  experimental 
spectra.   It  was  possible  to  obtain  the  frequency  of  rotation  as  a 
function  of  temperature  from  this  method,  and  an  Arrhenius  plot  enables 
the  determination  of  the  AH  and  AS   for  internal  rotation  of  the  methyl 
group.   The  results  are: 

AH  -  3.0  +  0.3  kcal/mole 

AS   =  6.0  +  0.3  cal/mole  deg 
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The  predicted  and  experimental  spectra  for  a  constant  magnetic  field 
orientation  and  varying  temperature  are  shown  in  Fig.  11A.   The  dotted 
line  spectra  are  the  calculated  spectra.   Figure  11B  is  the  Arrhenius 
plot. 

C.  LOW  TEMPERATURE  OBSERVATION 

Spectra  taken  for  single  crystals  of  strontium  acetate  hemihydrate 
irradiated  and  observed  at  liquid  nitrogen  temperature  without  warmup 

show  a  multiline  pattern,  as  can  be  seen  in  Fig.  12,  which  does  not 

o 
show  up  again  after  the  crystal  has  been  warmed  to  -100  C  and  then  re- 
cooled  to  -160  C.   The  orientation  of  the  magnetic  field  is  <f>  =   45  , 

e  =  70°. 

D.  MOLECULAR  ORBITAL  CALCULATIONS 

Molecular  orbital  calculations  for  CH  CO   CH  CO,  and  CH  CHO"  based 
on  the  CNDO  approximation  were  made  [27,28].   For  CH  CO  ,  the  calcula- 
tion predicts  the  unpaired  spin  density  to  be  localized  in  the  two 
oxygen  2p   orbitals,  with  0.499  spin  density  in  each.   For  CH_C0,  the 

£t  -J 

calculation  predicts  an  average  spin  density  on  the  methyl  carbons  of 
0.01916  which,  when  multiplied  by  500  gauss  which  is  the  hydrogen  atom 
coupling  constant,  yields  a  9.6  gauss  coupling  constant.   The  spin  den- 
sities are,  for  the  carboxyl  carbon,  0.57;  on  the  oxygen,  0.32  and  on 
the  methyl  carbon,  0.09.   The  calculation  for  CH  CHO  reveals  that  the 
average  spin  density  on  the  methyl  hydrogens  is  0.014076  which  gives  a 
coupling  constant  of  23.4  gauss.   The  spin  density  on  the  oxygen  is  0.35 
and  on  the  carboxyl  carbon,  0.47.   The  isotropic  coupling  constant 
calculated  from  McConnell's  relation  [l]  for  the  inequivalent  hydrogen 
is  10.7  gauss. 
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FIGURE  11  B.   Arrhenius  Plot  for  determing  activation  Energy 
for  rotating  methyl  group  in  -100°  C  species. 
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E.   SUMMARY  AND  CONCLUSIONS 

Radiation  damage  of  strontium  acetate  hemihydrate  single  crystals 
initiates  a  multistep  solid  state  free  radical  chemical  reaction.   At 

liquid  nitrogen  temperature  to  -100  C,  a  radical  of  the  form  CH   -  X 

CH 
with  a   3  =  9.2  gauss  and  a  pattern  of  eight  equally  intense,  equally 

spaced  lines  exists.   The  lines  exhibit  exchange  effects  to  become  at 

-100  C  a  four  line  1:3:3:1  pattern.   The  exchange  is  attributed  to  a 

rotating  methyl  group  with  barrier  to  internal  rotation  equal  to  3.0 

kcal/mole.   At  room  temperature,  the  species  is  of  the  form  CH  CH-X 

CH  H 

with  a   3  =  about  25.2  gauss  and  a  =  about  20.  gauss.   The  spectrum 

of  the  room  temperature  species  is  two  1:3:3:1  patterns  offset  from 
one  another  because  of  splitting  by  a  hydrogen  atom.   In  addition,  a 
broad  line  species  exists  at  room  temperature  which  may  be  a  triplet 
or  a  strontium  ion. 

Comparison  of  the  above  results  with  the  molecular  orbital  calcula- 
tions shows  that  the  methyl  group  coupling  constant  of  9.6  gauss  calcu- 
lated for  CH  CO  compares  favorably  with  the  9.2  gauss  methyl  group 
coupling  constant  determined  experimentally  for  the  low  temperature 
species.   CH„C0  would  also  be  expected  to  give  the  same  spectrum  as 
the  low  temperature  species.   The  calculated  methyl  group  coupling  con- 
stant of  23.4  gauss  for  CH„CH0  compares  favorably  with  the  approximate 
25.2  gauss  coupling  constant  for  the  room  temperature  species.   CH„CH0 
would  give  the  pattern  described  for  the  room  temperature  species. 
CH  CO-  is  ruled  out  as  a  possibility  because  of  its  spin  density  distri- 
bution on  the  two  oxygens.   Confirmation  of  the  identity  of  the  species 

in  the  crystal  might  be  made  by  irradiating  a  single  crystal  of 

13 

Sr(CH        CO    )    *%H  0   and    subsequently    examining    the  ESR    spectrum. 
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APPENDIX  A 

A  MATHEMATICAL  MODEL  FOR  THE  MOTIONAL 
AND  EXCHANGE  NARROWING  OF  SPECTRAL  LINES 


The  temperature  dependent  behavior  of  the  spectra  of  free  radicals 
with  a  rotating  methyl  or  -CH  group  may  be  predicted  by  a  computer 
program  based  on  a  mathematical  model  developed  by  Anderson  [ll]  and 
modified  by  Sack  [12]. 

The  phenomenon  of  narrowing  of  magnetic  resonance  line  breadths 
takes  two  forms;  motional  narrowing  and  exchange  narrowing.   The  mecha- 
nism is  in  principle  the  same  in  both  motional  and  exchange  narrowing. 

The  model  herein  developed  assumes  that  the  precessing  moments 
give  rise  to  a  radiated  electromagnetic  wave  which  is  undergoing 
frequency  modulation,  because  the  magnetic  interactions  act  to  change 
the  frequency  of  the  precession.   The  frequency  modulation  is  chang- 
ing in  a  random  way  in  time  due  to  the  effect  of  the  non-magnetic 
motions  on  the  magnetic  interactions. 

It  can  be  shown  in  general  that  the  spectrum  of  the  radiation  from 
any  quantum-mechanical  system  is  given  by 

°°  -iiut        2 

(1)  I(u>)    =  Trace    |\      M-(t)e  dt|  where  u.(t)    is    the   radiating 

—  03 

dipole  moment  matrix  in  the  equation: 

(2)  if)  -r^  =  Hu.  -  |i,H  where  H  is  the  Hamiltonian. 

In  any  problem  involving  the  type  of  narrowing  dealt  with  here  the 
Hamiltonian  may  be  split  up  into  three  parts: 

H  =  H  +  H  +  H 
o     p    m 
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These  three  parts  are  the  unperturbed  Hamiltonian,  the  perturbing 
Hamiltonian--generally  just  the  dipolar  interactions  between  the  moments- 
which  does  not  commute  with  H  and  thus  can  change  the  frequencies  radi- 
ated by  the  system  over  some,  known  range,  and  the.  motional  Hamiltonian 

whose  characteristic  is  that  it  commutes  with  both  H  and  u, ,  and  thus 

o 

can  have  no  direct  effect  upon  the  radiation  emitted  or  absorbed  by  the 
system,  but  does  not  commute  with  H  and  thus,  by  the  relation 

(3)  ihE      =  [H,H  ]  =  [H  ,H  ]  +  [H  ,H  ] 

p       P      o  p      m  p 

brackets  denote  the  commutator  can  cause  a  time  dependence  of  H  .   It 

P 

is  this  time  dependence  which  narrows  out  the  line  broadening  which  H 

P 

otherwise  would  cause. 

The  random  frequency  modulation  picture  of  the  narrowing  process 
can  be  derived  from  the  following  physical  assumptions: 

(a)  equation  (1) 

(b)  equation  (3) 

(c)  H  |JL  -  U.H   =0 

m      m 

(d)  H  H   -  H  H   =0 

mo    o  m 

(e)  and,  in  addition,  the  assumption  that  H   is  small  enough  that 

P 

it  has  not  important  matrix  elements  connecting  different  states  E. 

and  E  .     of  the  unperturbed  energy  H  . 

In  any  real  substance,  assumption  (e)  will  not  be  rigorously  true,  but 

may  be  looked  upon  as  an  approximation  which,  if  carefully  handled, 

leads  to  good  quantitative  results.   Under  the  above  assumptions,  the 

spectral  intensity  (1)  will  be  computed.   Consider  one  spectral  line  so 

that  only  the  element  of  U.(t)  which  connects  two  unperturbed  levels  of 

H  E.  and  E.  need  be  considered.   Then 
°   i      J 
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(4)         I(u>)  =  |C   ^..(t)e"iUJtdt|2 

i-j      J-=o   LJ 

Since  that  part  of  the  time  dependence  of  p, .  .  (t)  due  to  the  unperturbed 
Hamiltonian  is  known,  the  following  transformation  may  be  performed: 

i         ,       S  /•        (0)   v 

P».  .  =  P.  .(t)  exp(iu>.  .  v  yt) 


where 


Now 


(0)  =  E.  -  E. 
ij 


ift  -^  =  [H  +  H  ,  p]  =  H  p'  -  p'  H 
at     m    p        p         i 


because  of  assumption  (a)  and  the  fact  that  the  transformation  has  re- 
moved the  terms  of  the  time  equation  (2)  for  p  which  depend  on  H  .  By 
use  of  assumption  (e)  ,  p .  .(t)  may  be  calculated  as  follows: 


where 


Then 


ift[i.!  =  [H  p.1  -  p/  H J,, 
ij     P         P  XJ 


=  (H  ).  .p!  .  -  p.!  .  (H  )  .  . 
P  ii  ij     ij    P  JJ 


p!  .  =  -  iAu).  .(t)p|  .(t) 
ij        ij     iJ 


H  (t)..  -  H  (t).. 
Au>..(t)  =   P   1X       -B--LL 


ij 


fi 


p.!.(t)  =^i.(°)  exp(-i^   Au>i.(t,)dt') 

This  is  inserted  into  the  intensity  formula  (4)  which  gives,  for  the 
lineshape 

.c00  (o)         rfc  i2 

(5)      I..(U>)~  \   exp(-i(a)  -  (i).  :  ;)t  -  i\  Aco.  .(t')dt') 


8i 


Ato   (t)  is  a  random  function  of  time;  its  value  at  any  time  depends  on 
ij 

the  values  of  the  diagonal  elements  of  H  at  that  time,  and  these  will 

P 

change  in  a  random  way  at  a  rate  controlled  by  the  motional  Hamiltonian, 

as  is  seen  from  the  time  equation  (3)  (there  is  no  time  dependence  due 

to  H  because  these  are  diagonal  elements)  : 
o 

iftH  =  [H  ,  H  1 
P     m    p 

The  reason  for  assuming  randomness  is  that  in  the  important  cases 

the  effect  of  H  back  on  H  can  be  neglected,  so  that  the  motions 
p  m 

embodied  in  H  appear  to  the  magnetic  quantities  to  be  uncorrelated. 
m 

Equation  (5)  and  the  idea  of  Act),  .(t)  as  a  random  function  of  time  are 
n  i  j 

the  basic  ideas  of  the  random  frequency  modulation  model  of  narrowing 
phenomena.   The  term  Aco..(t)  is  simply  a  construct,  a  single  function 
which  is  meant  to  embody  all  of  the  broadening  effects,  diagonal  as 
well  as  off  diagonal.   In  order  to  find  the  line  shipe,  the  assumption 
is  made  that  Aco..(t)  can  be  treated  as  a  Markoffian  random  function. 

Equation  (5)  may  be  simplified  in  the  following  way:  Redefine 
co  as  to  =  co  -  co.  .    '  making  co  the  frequency  deviation  from  the  line 
center.   Introduce  the  phase  deviation  T|(t),  which  is  the  integral  of 
the  frequency  deviation: 


;0 
Then  (5)  becomes 


Tl(t)  =  9  Aco.  (t')dt' 


|\"  e"(lUt  -  "(t))dt 


Because  of  convenience,  the  correlation  function  form  of  the  Fourier 
integral  will  be  used.   It  can  be  shown  that 
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x/    \         0        iuuTj   0        -i(Tl(t'+T)    -    Tl(t'))J    , 
I(uo)    =  \      e        dT\      e  v  J  v      /ydt' 

or, 

00 

(6)  I(iu)    =  \     e±m  cp(T)dr 

^-00 

9(t)  is  called  the  correlation  (or  autocorrelation)  function  and  can 
be  written  as 


,  N    ,  i(Tl(t-h-)-Tl(t)v 
(7)  cp(t)  =  (e  y)      ave.  over  t 


=  <exp(i\    Aco.  .  (t ')dt '))  ave  over  t 


This  represents  the  averaged  "memory"  of  the  function  (j, !  .  at  the  time 

t  +  T  for  what  its  state  was  at  a  time  T  earlier.   The  problem  is  to 
find  cd(t)  under  the  assumption  that  Au>..(t)  is  Markoffian. 

A  Markoffian  function  is  a  function  for  which  the  probability  of 
a  given  value  f   at  the  time  t,  if  the  value  was  f     at  t  -  At,  is  inde- 
pendent of  the  value  of  the  function  at  any  earlier  time  than  t  -  A  t. 
Thus  the  probability  depends  only  on  the  value  of  the  function  at  the 
earlier  time  and  not  on  its  past  history.   The  notation  used  is 

Probability  =W(f  If  ,At) 

If  the  function  is  not  to  vary  in  a  senselessly  rapid  fashion  it  is 
known  that  W  must  be  proportional  to  At  for  small  At  unless  f   =  f  . 
Therefore,  W  may  be  written  as 


(8)        W(/1|.f2>At)  =  l(fvf2)   +   TT</1,f2)  At 

where  rr  is  the  probability  per  unit  time  of  a  transition  from  f      to  f  . 
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To  solve  the  stated  problem,  a  matrix  method  of  solution  is  used 
in  which  the  probabilities  W  (f.,/-)  are  considered  as  two  dimensional 
matrices  connecting  the  sets  of  indices  /   to  f  .   The  assumption  is 
made  that  Ato  has  a  discrete  set  of  values. 

The  average  in  (7)  may  be  written  in  another  way.   The  interval  T 

T 

is  divided  into  n  equally  spaced  steps,  such  that  —  is  very  small  com- 
pared with  the  rate  of  change  of  Au>.   Then 


n       t  rt+T 

(9)  exp        iS       -  Au>(t  +  m^  )  =  exp        iV  Au)(t')dt 

m=l  J+- 


The  probability  that  A(o  takes  on  a  certain  set  of  values  at  the  various 
instants  in  (9)  is 

(10)   W1(Au)1)  W(Au)1|Aco2,J)xW(Aa)2|Aau3,J)x  •••  xfl(A»n-1|  Au>n>  J  ) 

cp(T)  may  be  found  by  summing  the  product  of  (9)  and  (10)  over  all  the' 
possible  combinations  of  Acu's  since  these  are  respectively  the  value 
and  the  probability  of  a  given  averagand 

n  n 

.     T 


cd(t)    =  S    S     S     •••    (exp(i  -^  E     Ao>   ))   xW.  (Auk  )   it     W(A(u      .to,-) 
a...    a...    a...  _t      ni   '  _„  m-  L        m     n 

(ID 


Ao)  Aci)  Au)  m=l  '  m=2 


n-1 

=  S    S  •  •  •   S  tt      (exp(i-  Auo  )   x  W(Auu    I  Auo    ,,,—))    exp    (iAuo  — ) 
.       A  ,  n        m  m '      m+1   n  nn 

Au)  Ao)  m=l 

The  central  product  has  the  form  of  a  matrix  product.   It  will  have  this 
form  exactly  if  the  diagonal  matrix  Atu  is  defined  as  follows: 


Aid.   =  6(Ao).  ,Aoj.)  x  Ato. 
=ij       i   J      i 


The  probability  matrix  is 


W  (-)..  =W(A(W.|Aa).,  ") 
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W  (Aw)  is  a  row  vector  in  Aw  space,  and  _L  (Aw),  the  vector  all  of  whose 
components  are  unity  is  a  column  vector. 

T  T 

Expanding   exp(iAw     — )    and   neglecting   Aw     —     with   respect    to   unity, 

rnn        °      °nn  r 

(11)  is  written  as  a  matrix  product  multiplied  on  left  and  right  by 
row  and  column  vectors  respectively: 

(12)  cp(T)  =  W   •  (exp(i  -   &o)  ^(2)  )n_1  •  JL 

T 

Since  —  Aw  can  be  made  as  small  as  is  necessary,  it  is  permissible  to 
let 

T  T 

exp(i  —  Aw)  =  1  +  i  —  Aw 

n       —     n  = 

while  from  (7)  the  matrix  W  is  rewritten  in  terms  of  matrix  rr 

S<5)  =  i  +  eI 

This   gives    for    (12) 

(exp(i  I  Aw)   W(^))n_1   =   (1  +1    (i^u)  +n)    )n_1 


which  if  n  goes  to  infinity  is 

exp(T(i  Aw  +  rr)  ) 
so  that  (11)  becomes 

cp(T)  =  W   •  exp(r(iAw  +  rr)  )  •  I 

or 

cp(T)   =  S    W.  (Aw.)   iexp(T(iAw  +  tt) V   Aw  Aw 
i,j  _i       i      L  =       ='J         i      J 

(13)  =  W  expO(iQ  +  ^  )]    •    I 

where  the  n  components  of  the  vector  W  are  proportional  to  the 
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occupation  probabilities  of  the  sites  in  equilibrium,  1  is  a  vector 
with  all  components  equal  to  unity,  Q  is  a  diagonal  matrix  with  ele- 
ments cu .  and  rr  is  a  matrix  with  elements 
J 


TT.,   =  p 


TT.  .  =  -S  P., 


j  ^  k    p  =  probability 


For  cp(T)  given  in  (13),  equation  (6)  may  be  integrated  as  follows: 

\    exp(-iU)T)cp(T)dT   = 
Jo 

T 

W    •    [expiT(TT  +  m  -    iu>E   )>]      (tt  +  iO  -    icoE    )"      -1 
—  L     =  =  n   jo—         =  n  — 

where  E   is  the  unit  matrix  of  order  n.   For  T-»<»  the  integral  converges 
only  if  the  real  parts  of  all  the  eigenvalues  of  i  Q  +  tt  are  non- 
positive.   From  physical  arguments  it  is  clear  that  the  integral  must 
converge,  so  that 

(14)     I(u>)  ~  2  Real  W  •  (iuaE   -  iQ  -  tt)"1  •  1 

The  calculation  as  to  how  the  shape  of  a  composite  spectral  line 
depends  on  the  transition  probabilities  between  the  sites  is  thus 
reduced  to  the  inversion  of  a  matrix  with  complex  elements.   A  computer 
program  to  perform  the  calculation  in  (14)  has  been  written  and  begins 
on  the  next  page. 
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THIS  PROGRAM  CALMJLATEC 
rpMFpcITc  SPFCTRAL 
RFFForNjCF0    FOP    THF 


THF  INTFMSITY  DISTRIBUTION  OF 
LINF  WTTH  PXCHANHE  NARROWING. 
THEORY    APC 


J    PHYS     ^HC     JAPAN, °,88P 
M^LPfMl.40     PHYSICS ,1 ,1*3 


WRITTFN    PY    PHRFPT    G.     GTSCH 


ADPTL        1969 


1*T   CAPO       sft    1    PnR    p^P  ?    POR    NMR/NO.    OF    PLOT*     pfp 

GPAPH/Nri.     OF    OP&PHS/NO.     OF    POINTS    TO     PE    PLOTTED/ 
WIDTH    PC    FRFOUENCY    SCAL^    IN     I NCHF~ / DI "ENS I nM    OF 
m/trtx/NT.     HF     NON7FRO    OFFfJIAGONAL     ELCMENTS     ABOVE 
njAGrw L/LI MF    WIDTH    TM    INCHESITHE    RECIPROCAL    OF 
T2)/THF    HFir.HTfIN    INCHES, OP    THF    HIGHEST    ABSOoeTipM 
RF*-K  11,12,12, 13, 312, 2FIO. 5 


2ND    fipD 


VFCTOP     ^F     INTFNSTTIF? 


6F10 


5 


3on 


CAPO 
LINES 

INCHF 


ptjgonal    matrtx   m=    frpqiif^CIf^    rp    p=*pnancf 
MATRIX    ELFVCNT    *NO    POSITION    nc    LINE    IN 
6 (211, FIG. 5) 


4TM    fAPH  THE     TRAMSITTTN    PROBABILITIES    TR£N"FlPMCp    JP 

'HE    IMFHFS    CC*LF    *c    IS    thf    pthfq    hata.       tujE    first 
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TCI  I 
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CAPO  1  5A  4 

7TH    fAPD  ADDITIONAL    INFORMATION.        TY'F     IN'     Ac     MANY 

Ci*PDc    t*    THFPP     APT    GRAPHS  **° 


IMPLICIT    CC*FLFX*16(U,0,V,X> 

IMPLICIT    pcAL*p( 4-H,n,P, b-t, 
RFAl*^     XPEO,XINTFN,YT NTFM 
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•  / 
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1  •  / 
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nn     icn 1     JT =1 fNMfo 
on    1001     T  1^1 f | i  l 

c=t    up    ?prBABILITY    MATRIX 

nn    cr     I  =.]  ,  mm, 2 
T"l=T-i-l 

I  *  =  T I  R  (  T  ) 
IP=T/R(  t  i  ) 

rr    "PPR(  T '  ,!P)  =-RFRpp  ((JT-1)*!IL+TJ) 

nn    *c    t  =i f m 

n°   5«   j=1,n 

cc     ppnp(    |  f j )=oorp  (  I  ,J) 

O1"*     22     1  =  1  ,N 
2?    oppQ (  T  f I  )=C.C 
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UPP.nPU  t  J)=UPPDR(I  f  Jf*(-1.0v0«0) 

'c  uu(  t  , j)=uo*cg*  ( i  ,;n+npon3(  i ,  n 

nn   loop    wf-TFH   CALCUL/STFS    IMT  =  NSITTC*^ 
nn    lOOC    M=itNUMPTS 

no   li    1  =  1  vn 

0"    11     J=1,N 

11    U( I ,J)=UU( T,J) 

|)prcp(  ^)-ppcci(  M) 

UFPPn( M)=UFRFQ(M) *<C.C,1.0) 
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PLFMFNTS 

OP    2    T=1,N 
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<:=o 
INTENSITY    TS    TWKF    THF    ORL    PART    CF    THIS    PFS'JLT 
1000    7  INTFMf  *')=2.  *c 

FisjD   C&LCUL^TIONS.       STAPT    cnMPijTj^p,   ppa,phc 
IF( KKL.FQ.l)    GO    to   450 

DO    5CC    !=1,M 
XOFOf  I  )=ppco(I ) 
5C0    X!*'T(f  N(  T  )=ZINTFN(I  ) 

Icf  J  I  .GT.  1  )    f,n    rr.    cqn 

ipf  tj.pt.  i )   r,n   th   oqp 

Tnp=r.o 

r><~    4^    T=ltNLMCTf 

TF(  TPP.GT.XTNTFN(I  )  )     G^    Tp    444 

TriP  =  XINTFN( T ) 
4^4    rOMTTMUF 
qcq     pn     4C4     I=ltN(JMPTC 
6  5^     XTMTFNf  I  )=XT  NTPM(  I  )  *?,rA  Lp/Tr,° 

1  F  f  IJ.FO. 1)00    TP    765 
On     76f     T=7T12 

7f6     ITL( T )=PL*NK(I ) 

on   tp   7^7 

765    ^n     7AP     l=7tl 2 

7A9    TTL(  T  >=ITITLF<JI  ,1 ) 

167    CONTINLF 

LB (  =1.  &RFL(  (  JI-1  )*LILM  J) 

TFf LLL.FQ.l)     OP    Tp    204 

tf<  Tj.Fn.  i)    pn    tp    201 

Ic( IJ.rO. ILL)    00    TP    202 

CALL     OP£W(N"J"DTCtxiNTCNfXPPQ,?  ,0,L8L,ITL,1.C,  1.0,0,0, 
12,2,c,THinH,0,UAST) 

OH    Tn     1001 

201  P£LL    DRAW(V'UMPTS,xiNTc:NtxPFO,l,0,LpL,ITL,l.C,  1.0,  0,o, 
12,2,°,IHlOHt0,L AST) 

on    tp    icoi 
20*    TALI.    np^wfNll^P1"^,  XINTFNfXo.PO,OtO,LPLtITL,  1.0,  1.0,0,  C, 
12,2,9,IHIGH,0,LAST) 
on   tp    icri 

202  PAI.L    DP/W(NUMPTC, XTNTPM,XP^O, 3, o,LRL,ITL, 1.0,  1.0,0,0, 

12,2,9»IHIGHf0,L*ST) 

On    tp    iCCi 

FS"     SFCTJPN 
450    on    P0C    1  =  1, nvi 
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CUCDUTFC     NMIMPRICAL     DERIVATIVE 


1C* 
5  c  1 


3f  6 

SA7 


301 


30  < 


30? 

10C1 
100 
102 
105 
lin 

120 
121 


HIMTF 

r-o    Rt 

X«?rO( 

YTMTp 

T  F  (  J  T 
IF<  I  J 

T0P=C 
n^   44 

I*-(Tn 

Tip-;, 
rrMJT 
nn    ^c 

YTNJTT 
TFf  T.I 

OP    86 

TTL(  T 

C,n  to 
no    of. 

TTL(  T 
CHN'TI 
LBI.=L 
Ic(  LI 
I F  (  I J 
TF(  IJ 

cai  l 

1  ^  t  I  H  T 

GO    t 
CALL 
I3f THT 

C  6  L  L 

1  <? , !  H  T 

GO    TP 

C/l  L 

ICtlHI 

cpmtt 

pno  m^ 

FORMA 

FHRMyv 

FP.PMA 
FORM/ 

FOPM,». 


N(  T  )  =  (  7INTFM(  1  +  1  )-7TNTCM(T  )  )/  7pJC 

1     T=1,M 

T  )=focQ(I  ) 

Ml  T )=HT  NTFNI I  ) 

•  GT.u   r,n   Tn   Q9C 

•GT.l)    GO    TO    QQO 
.C 

".fJ.A^I  YTNTFNIT  )  )  )     Pn    TP    44* 

RSI YlNTrN( I  )  ) 

M|  IF 

c      *=1,NM1 

N{  t  )=YT\TFMU  )  *cr*LF/  (TOP*?.  H-4.5 

.FO.IJGO    TH    86K 

*  T-7,12 
)=BLANK(I) 

Rf  7 

Q    T=7tl2 

)=TTITLF(  JI  ,1  ) 

M  IT 

*RrL(  (.11-1  )*LIL+I.J) 

L.F0.1)    GO    Tr    30^ 

.FO.l)    n,o    JO    301 

.FQ.LLL)    GO    to    3Q2 

DPAV(NMl,YINTFN,XRFQ,2,Q,LBL,ITlfl .0 , 1.0,0,0, 2, 2, 

GH.O.U  rr) 

IC01 
DRAW(NM1 ,YINTEN,XRFO,1,0 ,LRL, ITL, 1  .C , 1.0,0,0, 2, 2, 
GH,C,LA$T) 

1GC1 
DOAW(NMltYTWTFN,XDFOf0,0,LRL, ITL, 1.0,1.0, 0,0, 2, 2, 
OH ,0, LA  ST) 

1001 
OPAMNMi  ,YINTFN,XRFO,3,0,LRL,  ITL,  1  .0  ,  1.0,0,0,  2,  2, 
GH.CLASTJ 
NUF 

T(  II, 12, 12, 13, 312, 2F10. 51 
T(6( 2I1,F10.51I 
T(7C(  ID) 
T(6F10.5) 
T(  le£4) 
T(6£8) 


11 
12 
13 


14 


SUP  TUT  IMF  GAUSS4(N.FPfA  fX,KER) 
IMPLICIT  CPMFLFX*1<S  (  A-0,F-HtP-Y) 
niMPMSIPN    M20, 20)  ,X(20,?0) 


OP     1     I=1,M 

nP     1     J=1,M 

1 

X(  I  ,J)=(0. 0,0.0) 

Dr     7    K=ltN 

2 

X(K?K)  =  ( 1.0,0.0) 

10 

nn    3A    L=1,N 

KP  =  C 

7=0.0 

nn    12    K=L,N 

IF! 7-CDAB«(A(K,L))  ) 11  ,12,12 
Z*COABS(A{K,L)l 

CHNTINUF 

IF( L-KP) 13,20,20 

nn    14    J=L,N 

S=A(Lf J) 

A  (  L , J ) =A ( K© , J ) 

A  (  K  P  ,  J )  ■  S 

nn    i*    j=i,n 

S=X(L,J) 

X(L,J)=X(K»f J) 
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15 

XCKP, J)= 

?0 

30 

!F( L-N)3 

31 

LP1-=L  +  1 

no   3*   k= 

T  F  (  C  P  *  B  c 

32 

r; tip=£ ( 

on    33    .1  = 

33 

M K, J)=* 

nn    ^s    j  = 

3^ 

X(  K,  J)  =  X 

36 

rnMTjrMdF 

34 

mMTiM.'P 

*C 

nn    4  3    1= 

II=m+1-I 

on  43    J= 

*  =  ( o.o,c 

I F  (  1 1  -  N  ) 

M 

I  TPl  =  !T  + 

nn   4?   k  = 

42 

c  =  c+^ (  II 

A3 

Xf  !!  ,J)  = 

k-FD  =  l 

RFTljOM 

50 

KFQ-2 

RFTHPN 

CMP 

IML.LI  )-F°l*0, 50,30 
1 ,34,34 

LP1,N 

(MK,U)  )32,36,32 

K ,  L )  /A  (  L  ,  L ) 

L  p  1  ,  N 
(KtJ)-PiTTn*A (L, J) 

ItN 

f  Kf,l)-P4Tin*x(L,  J) 


i,n 

1,N 

.0) 

41,43,' 3 

1 

1 1  P 1 ,  N 

.K)*X(K,J) 

( X( 1 1 » Jl-5) /*  C I  I  , T  T  > 
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